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^ I Abstract. 

We investigate some implications of the freezing scenario proposed by Carpentier 
and Le Doussal (CLD) for a Random Energy Model (REM) with logarithmically 
correlated random potential. We introduce a particular (circular) variant of the 
model, and show that the integer moments of the partition function in the high- 
temperature phase are given by the well-known Dyson Coulomb gas integrals. The 
CLD freezing scenario allows one to use those moments for extracting the distribution 
of the free energy in both high- and low-temperature phases. In particular, it yields 
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( the full distribution of the minimal value in the potential sequence. This provides 

an explicit new class of extreme- value statistics for strongly correlated variables, 
, manifestly different from the standard Gumbel class. 
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1. Introduction. 

The Random Energy Model (REM) introduced by Derrida in [I] is characterized by the 
partition function Za = J2i=i e~^ Vi , where f3 is the inverse temperature, and Vi for i = 1, . . . , M 
are random variables with the typical variance < V- 2 >= V 2 = 0(log M) for M — > oo. Such a 
model, as well as its generalized version (GREM) [2j [3] continues to play a paradigmatic role in 
statistical mechanics of disordered systems. Simple enough to allow for a detailed analytical 
investigation by various methods, the freezing transition exemplified by REM appears to be a 
rather generic phenomenon. It emerges with surprising regularity in a variety of physical 
situations, ranging from transparency of random media [I], directed polymers in random 
environment [5], p-spin glass models and the glass transition [6], random heteropolymers and 
models of protein folding \1\ to properties of quantum particles in a random magnetic field 
[HJ [9] , and thermodynamics of a single particle in random Gaussian landscapes [10] . It is also 
a rich model from purely probabilistic point of view [11], and has an interesting dynamical 
counterpart: aging [12]. 

As the low-temperature behaviour in statistical mechanics is obviously controlled by 
the lowest available energies, it is not surprising that a detailed description of the freezing 
phenomenon is intimately related to the so-called extreme- value statistics of random variables 
[13j . For independent, identically distributed variables Vi, the cumulative probability 
distribution P m (x) of V m = min{l/i, . . . , Vm} is well-understood in mathematical literature 
and, provided the support of the distribution extends to — oo but decays faster than any 
power-law, is given in the limit M > 1 by the Gumbel law: 



Here the constants om and bu depend explicitly on the distribution of Vi, but the double 
exponential shape of the Gumbel law is very robust (universal). This universality extends to 
a very broad class of correlated variables provided those correlations decay fast enough, see a 
detailed description in p3]. For a mathematically rigorous analysis of REM and GREM based, 
in particular, on the extreme value statistics, see |15| ; see also [TU] for recent developments. 

A few years ago Carpentier and Le Doussal (CLD) [9j studied a specific case of correlated 
random variables that is arguably the richest, most challenging, and relevant for applications. 
More precisely, thinking of the index i as referring to the sites of a certain lattice, the covariance 
{ViVj) considered by CLD depended on the distance d{i,j) between those sites logarithmically. 
An independent support of the fact that logarithmically correlated potentials should play a 
special role was obtained recently through a thorough analysis of statistical mechanics of a 
single particle in high-dimensional random energy landscapes |10| . 

To understand the extreme value statistics in the logarithmic case and to relate it to 
a REM-like freezing transition CLD developed a powerful, albeit non rigorous, real space 
renormalization group approach to the distribution V(Z) of the partition function Zq. In this 
way they discovered that logarithmic models represent, in a sense, continuous analogues of 
directed polymers on disordered trees, a result somewhat anticipated in [SJ. The statistical 
mechanics of the directed polymer problem on a tree is known to be amenable to a travelling 
wave analysis, see the celebrated paper [5]. Similarly, for the logarithmic case the Laplace 
transform of the distribution V{Z) was also shown to satisfy a kind of travelling wave 
equation (see also works [T7] for a relation between the travelling waves and the extreme 
value statistics, in particular, in the context of the zero-temperature directed polymer problem 
on a tree). Solutions of equations of that type are known to exhibit a characteristic change 
of the shape at some critical front velocity, and following [5j that change was interpreted 
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by CLD as a signal of a REM-like freezing. In particular, the CLD analysis revealed that 
such a transition implies, among other properties, a universal non-Gumbel shape of the far- 
left tail for the cumulative distribution of the minimal value of logarithmically correlated 
variables: P m (x — > —oo) ~ 1 — Const \x\ e ax which is clearly different from the Gumbel tail 
Pm{x — * — oo) 1 — Const e ax . 

Although CLD's renormalization group is able to predict successfully the universal far-tail 
features of the distribution expected to be shared by all logarithmically correlated potentials, 
the calculation of the full distribution for a given potential is beyond the scope of that method. 
Indeed, the renormalization procedure needs specifying what the authors called a "fusion of 
environments" rule [9]. As CLD convincingly argue, the precise form of that rule is not 
important for recovering universal properties. But the actual shape of the travelling wave 
equation certainly depends on the particular fusion rule employed. As no guiding principle for 
the fusion rule selection had been provided, a more detailed analysis of specific models appears 
to be problematic. 

The present paper grew out of attempts to overcome the above difficulty. Our main 
observation is that for a particular variant of the logarithmically-correlated REM a more 
detailed analysis seems to be possible. In particular, we are able to conjecture the full 
distribution function P m (x) pertinent for the case considered. Namely, after appropriate 
rescaling the cumulative distribution P m (x) of the minimum turns out to be given by the 
central result of this work: 



where K\{x) stands for the modified Bessel (Macdonald) function, and the parameter (3 C 
is given by the inverse transition temperature in the model. As such, the corresponding 
distribution is manifestly different from the standard Gumbel double-exponential form. 

Our method draws its inspiration from CLD analysis and is essentially based on the 
pattern of the freezing transition revealed in [3 12]- However, we do not resort to the 
renormalization group construction or travelling wave technique, and in this way circumvent 
the need to know the microscopic fusion rule. Instead, we demonstrate below that assuming 
CLD freezing scenario provides a way to extend the moments of the partition function from the 
high-temperature phase to the region below the transition. Roughly speaking, if those moments 
can be explicitly calculated above the transition,- as is the case for both the standard REM as 
well as our variant of the logarithmic model,- they can be used to conjecture the shape of the 
Laplace transform of the probability density above the transition point. Then CLD approach 
can be used to recover the full distribution of the partition function/free energy below the 
transition, yielding in particular the extreme value statistics. 

The structure of the paper is the following. We start with a short general discussion of 
CLD freezing scenario, in particular its implications for the moments of the partition function 
in the low-temperature phase. As an illustration we show how CLD scenario can be used in 
the standard REM case to recover the exact low-temperature expressions for the partition 
function moments, obtained long ago by Gardner and Derrida via a rather tedious analysis [3]. 
After that we introduce and analyze the particular (circular) version of the one-dimensional 
logarithmically correlated REM. 

2. CLD freezing scenario: general relations and implications 

The central object of the subsequent analysis is the Laplace-transform Gp{p) of the probability 
distribution V(Zp) of the partition function: Gp(p) = (exp{— pZp}). Here and henceforth the 
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angular brackets stand for the expectation with respect the distribution of random variables 
V{. Our approach is based on the assumption that such a Laplace transform can be efficiently 
found in the high-temperature phase. 

Following [9] we introduce the variable x via p = e@ x and consider the function 
G(x) = (exp{— eP x Zp}). Extending the Derrida-Spohn scenario, CLD postulate that in the 
thermodynamic limit M»l the latter function has a shape of a travelling wave, that is 

G{x) =gp(x + mp(L)), mp(L > 1) « c(p)L + l.o.t., L = lnM (3) 

where we introduced the parameter L to identify with notations used in [9j, and l.o.t. stands 
for the lower order terms when L — > In general, both the travelling wave profile gp(y) 
and the wave velocity c(/3) depend on the inverse temperature (3. The REM-like transition in 
this approach is described by a "freezing" of both the velocity and the profile function at a 
certain transition temperature (5 = (3 C so that in the full low-temperature phase (3 > (5 C one 
has: 

9p(v) = SaM mp(M » 1) « c(f3 c )L + 1. o. t. (4) 

By the very definition of the function G(x) we then have the following relation for (3 > f3 c : 

(expi-f&Zf,}) = G{x)\p> Pc = g Pc (x + c(/? c )L) (5) 

which fixes the shape of the Laplace transform Gp>rj c (jp) = (exp{— pZg}) in the whole 
temperature range below the transition point. 

Equipped with such a scenario, we will exploit that the knowledge of the Laplace 
transform allows one to calculate the moments of the partition function below the transition 
by employing the standard identity: 

1 r°° 1 r°° /I \ 

W = ml dpp " 1{e ) = W)L ^"M/J lnp+C(ft)i )' (6) 

where T(x) stands for the Euler Gamma-function, and we used that x = 4 Inp. 

According to the freezing scenario, the mean value of the free energy in the low- 
temperature phase is to leading order in L ^> 1 temperature-independent, and is given by 
(F) = —jf(\n.Zp) = —c(/3 c )L. Being interested in the fluctuations of the free energy, we 
introduce the random variable f = F — (F) whose probability density we denote Vp(f). After 
changing the integration variable in © to y = jjlnp + c(f3 c )L, introducing s = —jiv and 
integrating once by parts we observe that takes the following form: 



f"0O 1 POO 

-*/\ = / p-sZ-p^fWf = / p- s v 



-i e- sf Vg(f)df = v e 

f J-oo ^ J > r(i-f) 1 



dy (7) 



Thus, the only function needed for investigating the low temperature phase is the shape 
of the travelling wave profile gp{y) at the critical point (5 = (3 C . For a general non-zero 
temperature C < (3 < oo one can extract the explicit form of the free-energy distribution 
Vp{f) by noticing that the analytical continuation s — > is converts the relation ([7]) to a Fourier 
transform which as we shall shortly see can be frequently inverted explicitly. A particular 
simple relation is obtained in the zero-temperature limit (3 — > oo where the free energy simply 
reduces to the minimum value of all random energies in the sample F — > V m = minj{Vi}. It is 
immediately clear from ([7]) that 

$mVp(f) = ~9p B tf) (8) 



| Establishing the precise form of those terms, both above and below the transition, was one of the 
central points of CLD analysis. It is however of no direct relevance for us in the present paper. 
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yielding a very general relation between the shape of the critical profile gp c (x) and the 
probability density V m (x) = ~^Pm( x ) of the fluctuations of the extreme values in the sample: 
'Pm(x) = -^gp c {x). 

Let us briefly demonstrate how this method works for the standard REM with i.i.d. 

cm ( 1+ 7&) lnM 1-1 
Gaussian sequence of v^. Introducing the notation Z^' = eV , the analysis of [9] 

demonstrated that everywhere in the high-temperature phase (3 < (3 C the Laplace transform 

is given by 

G /3</3c (p) = f°° e- pZ T(Z) dZ « e- pzm , < P Z<® <C 0(ln In M) (9) 

Identifying L = In M as in ([3]) we find from ([9]) and the correspondence G(p)\ 
gp(x + c(P)L) the REM travelling wave profile which is given by 

^(y) = exp{-e^}, c (0) = I + A . (10 ) 



p=e 



/3x 



c 



In particular, when approaching the transition point (3 = j3 c the shape and velocity of the 
travelling wave tends to the limiting values 

9pM=eM-e PcV }, c ^) = J c ( n ) 

According to the general discussion, the Laplace transform (e~ pZl3 ) = Gp{p) of the probability 
distribution of the partition function in the low-temperature phase (3 > (3 C can be found from 

G{x) = g Pc {x + c(P c )L) = exp{-e^ +2L } =}► G p>P M = e~ c ^ , C M = e 2L = M 2 . (12) 

where we again used the correspondence x = 4lnp. In turn, the last expression can be 
immediately employed to recover the (non-integer) moments of the partition function from the 
first of relations in ([6]). Substituting there Gpyp^p) from (fT2j) yields after a straightforward 
manipulation: 



-, r ( 1 + 



— v j 

(V> = W y r(i+% ' Za = e ^ lnM (13) 

valid as long as — f < The latter moments are, to leading approximation, precisely those 
obtained by Gardner and Derrida [3], and coincide, as expected from general arguments, with 
the moments of a totally asymmetric Levy stable distribution of index (3 C / f3 [13 [181 IS] ■ Finally, 
the above moments allow one to recover the distribution of the free energy fluctuations in the 
low-temperature phase of the REM, which seems not to be written explicitly in the literature. 
Namely, making an analytic continuation v — * is/ (3, and introducing / = — ^ In Z/Z a we 
notice that (|13|) takes a form of the Fourier transform of the probability density for /, see ([?])• 
Inverting that transform gives: 

r REM, f , 1 r r~"l 1 rfl I "W ■ d V (141 

v ^ u) -2iL e WTf) r V + ftj ds --5f£ — r (!-„!)■ (14) 

In particular, the zero-temperature limit (3 — > oo coincides with the general relation ([8l) between 
P m and 5^ c , as given by Eq. (|lip . which immediately yields the famous Gumbel distribution 
for the minimal energy. It is also evident that the far- left tail / — ► — oo of the probability 
density Vp EM (f) is of the Gumbel form everywhere in the low-temperature phase f3 > (3 C , 
again expected from general arguments [131 Q21 US] . 
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3. The circular logarithmic REM 

3.1. Definition of the model and the moments in the high-temperature phase 

Consider the lattice of M points positioned equidistantly at the circumference of a unit circle. 
Their angular coordinates are given by 6k = §y&, k = 1,2,..., M. With each point we 
associate a Gaussian random variable Vi, with position-independent variance (V 2 ) = V 2 and 
covariances chosen to be: 

-0i 



Cki = (VkVt) 



-g 2 In < 4 sin 2 



(15) 



For the consistency of the procedure we have to choose variance V 2 in a way ensuring 
positive definiteness of the full covariance matrix with entries V 2 5ki+(1— Ski)Cki- The condition 
amounts to V 2 > — A max , with \ max being the largest eigenvalue of the matrix C with entries 
Cki- The matrix C is by definition a circulant real symmetric, with zero diagonal. Hence, its 
eigenvalues are given by \ q = Yli=2 Ciiuj^ 1 , where uj q = exp{ 2 ^-} are roots of M — th degree 
from unity. The largest eigenvalue corresponds to q = 0, when uj q = 1. Then we have: 

-A/-1 



Ar. 



Using the identity: 



AI-l 
„2 \ " 



-g 



s 



^ In <M sin 2 — } = -a 2 In < TT 4 sin 2 — 



M 



n 



M-l 



sin 5f we see that A r 



-2g 2 InM, so we have to choose 



the variance to satisfy V 2 = 2g 2 InM + W, with an arbitrary positive W > 0. In what follows 
we assume that W = 0(1) when M 3> 1 and therefore can be safely neglected if we are 
interested only in the leading terms in the thermodynamic limit. 

We define the partition function for our model in the standard way through Zp = 
Yl,iLi£~^ Vi i with the goal to evaluate positive integer moments ([Z^]™) in the limit M S> 1. 
The expected value of the partition function is obviously independent of the covariance and 
is given by the standard REM expression, (Zg) = exp{lnM[l + (3 2 g 2 ]}, so the first nontrivial 
moment is 

M M 

{[Zp] 2 ) = ]T <exp{-2/3VU> + E <exp{"Wi + V l2 )}) 



e lnAf(l+4/3V) _|_ gin M2f3 2 g 2 



M 

E 



2 sin 



M 



n -2/3 2 9 2 



(16) 



Introducing 9k = irk/M for A; = 1,2,... 
Riemann sum, and that in the limit M 



, M we see that the second term in (|16p is a kind of 
■+ oo can be approximated as 



M 

E 

11^=12 L 



2 sin 



tt(*i 



^2) 



M 



-2/3 V 



M 2 

TT 



pIT — 7v/M 






2 sin^ 


Itt/M 





de 



e 



,2 InM r(l-2/3 2 g 2 ) 
C T 2 {l-I3 2 g 2 ) ' 

lnM(l+20 2 a 2 ) 1 

= 2 n 2 j 



7r 2/3^ 9 ^ 2 /32 g 2_ 1 ) 2 2/3^ 



2q2-1 ) 



2(3 2 g 2 < 1 
2/3 2 5 2 > 1 



(17) 



The first line in (|17|) corresponds to the convergent integral in the limit M — > oo which can be 
easily evaluated by reducing it to the standard Euler's integral of the first kind , see [22], p. 
898. In contrast, the second line is obtained by extracting the leading term of the divergent 
integral. Comparing now second "off-diagonal" term in (|16p with the first "diagonal" one, 
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we find that for 20 2 g 2 > 1 the diagonal and off-diagonal contributions are of the same order, 
whereas for 20 2 g 2 < 1 the off-diagonal contribution dominates. Finally, we arrive at 



([Zp] )|m»1 ~ < i nM(1+4/ 32 g2) 



e 21.M( WV )glzWl ; 2 g 2 < 1/2 



1 + 



7r 2^s^(2/32 g 2_ 1 ) 2 2^s^-l 



W > 1/2 <18) 



The case of a general positive integer moment can be treated along the same lines. 
Denoting Xi = e~^ Vi we use 

M n-1 , M-l M 

= E< + E^rniE E «-' + ••• as) 

h=l 1=1 * ' h=l «2=*i+l 

where in the second line we regrouped the terms according to partitions of the integer n into 
sum of nonnegative integers with length k (i.e. the number of nonzero parts) taking values 
k = 1, . . . , n. For example, partitions of the length k = 1 are sets {p±, . . . ,pm} with all but 
one pj equal to zero, and with the remaining nonzero integer taking the value p = n. The 
total contribution of those partitions is obviously Ya^=i x i L which is the first term in (|19p . 
Similarly, the contribution of all partitions of the length k = 2 is precisely the second term 
in the above expression, and so on. Finally, we perform the ensemble averaging of the above 
sum using the identity: 



k r /_/_• • \x -, ~20 2 g%l q 



Ik 

p<q 



2 sin 



7r(ij 



M 



(20) 



valid in the case of all different indices in the set i%, . . . , if.. In the limit M — > oo we then find 
by inspection the dominating terms. After manipulations generalizing those we performed 
earlier for n = 2 case we find the following general expression: 



<[^r>U/»l « | e lnM(l + nW) (l) / n> V p2 g 2 ( 21 ) 



where 



in(i3 2 9 2 ) = ^ r Mi r de 2 r ^ 3 . . . r de n n [2 S m % - g ))r 2/3V ■ ^) 

n n Jo h 1 Je 2 ££ q 

The explicit expression for the factor 0(1) in the second line of (|21|) is rather complicated, 
but is actually not needed for our purposes. Finally, using the symmetry of the integrand in 
([22]) and noticing that \e 2ie p — e 2t9q \ 2 = 4 sin 2 (0 p — 9 q ) one observes that ([22]) is a particular 
case of the so-called Morris integral (related to the famous Selberg integral) whose value was 
first conjectured by Dyson in his studies of the Coulomb gas problem (see a very informative 
historic account and further references in [19j): 

^ - ^ c * - ■ r * n ^ - ^ = • <*> 

The integral is clearly finite provided /3 2 <? 2 < 1/n < 1, and is divergent otherwise. The 
condition (3g < 1 defines the high-temperature phase of the model. Note a certain similarity 
between our calculations and those arising in the framework of the multifractal random walk 
model of Bacry, Muzy and Delour 
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The crucial point of our analysis is the ability to offer the explicit form of the probability 
density V(Z) of the partition function Zn = Z > which precisely reproduces the expressions 
for the moments ()21|) . It is given by 



nz) = { ; ' ;:;* , (24) 




where we defined Z* = e and introduced for (3g < 1 the two functions: 

^=w(l) exp -\{i) }' z ' = r(i-gV) ' (25) 

and 

>V J \/ 47r In M Pg Z J \2\nM) V ; 

To understand the structure of V(Z) notice that the growth rate of the moments in the second 
line of Eq. (]2ip dictates that the far tail of the distribution must be of a log-normal nature. 
This is exemplified by the choice (|26p for V>(Z). On the other hand, the first line in (|2ip and 
the expressions (|23p yield the probability density of the form V<(Z) in Eq. (|25p . The crossover 
value Z = Z* is determined from the requirement for the leading exponential terms in the two 
pieces of the probability density to match smoothly, i.e. V<{Z*) « V > (Z*) for M S> 1. The 
factor f(x) in (|26p is assumed to be of order of unity when its argument is of the order of unity, 
and is otherwise left unspecified. Finally, we verify in the Appendix that the choice of P(z) in 
Eq. ([2"3]) ensures the required change in moments (Z n )M^>i to occur precisely at n = l/g 2 (3 2 . 

At the next step we use the probability density (|24|) for evaluating the Laplace transform 
function Gp(p) in the high-temperature phase. A somewhat lengthy but straightforward 
calculation reveals that the log-normal tail V>(Z) gives for M 3> 1 a negligible relative 
contribution to the Laplace transform, as long as we keep finite the value pZ e < oo. Effectively, 
it means that for our goals we can assume the partition function Zp to be distributed with the 
probability density V<{Zp) given in Eq. (|25p . After a simple transformation of variables this 
implies a rather simple asymptotic formula: 

roo _ 

Gp(p) = (e-* z e)\ M>1 * i e-t-^^dt, a = (3 2 g 2 <l, (27) 

J o 

Using such an expression , we can, for example, easily calculate the mean logarithm of the 
partition function, hence the mean free energy: 



(In Z) = lim 







oo 

dpp e - 1 {e- pZ ^) 



lnZ e -ar'(l) 



so that the mean free energy is given by 

{F) = -h\nZ) = -(\ + (3g 2 )\nM-\\n\T(l-(3 2 g 2 )\ - (3g 2 T'(l), P 2 g 2 <l (28) 



The leading term yields the expected universal REM expression for the mean free energy valid 
in the high-temperature phase, the rest corresponds to system-specific corrections. Those 
corrections diverge logarithmically when approaching the critical temperature P = 1/g = @ c , 
signalling of the phase transition. Note, that the same result for the free energy can be 
recovered by the standard replica trick using moments (|2ip . 

The fluctuations of the free energy around its mean value can be easily recovered as well, 
using the explicit form of the distribution T , < (Zp). Namely, introducing 

/ = F- (F) = -±ln{Z/Z e } + /o, /o = -^T'(l) (29) 
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the equation (|25|) implies the following probability density in the high-temperature phase 

(3<(3 C : 

VpU) = ~ exp {£(/ - f ) - e^f-fA . (30) 
a [a J 



According to our previous discussion, a central role is played by G(x) = Gp (p = eP x 

Identifying L = InM we observe that G(x) = gp(x + mi) where mi = ^lnZ e . Using Eq. (f28j) 
we further see that 

m L \ L »i « l ~^L + 0(1) = c((3) L + O(l), c(/3) = ± + A , 

again in full agreement with CLD results in the high-temperature phase, with c{(3) interpreted 
as the travelling wave velocity, and the wavefront profile given by 



r°° [ e@ y 1 

9p(y) = J Q dtexpl-t- — \, a = 0Y<l. (31) 

3.2. Transition to the low-temperature phase in the circular logarithmic REM. 

To investigate the low-temperature phase for (3 > (3 C we rely upon the CLD freezing scenario. 
When approaching the transition point a = 1 / (3^ = 1 the profile (|3"T|) tends to a well-defined 
limit: 

g Pc (y) = 2e^ L K 1 (2e^ L ) , (32) 

where K\ (x) is the modified Bessel (Macdonald) function. According to the freezing arguments 
this shape via the relation (jHJ) is translated into the extreme value probability density ((2J) , which 
is our central result. This expression is non-Gumbel as the cumulative distribution behaves 
for x — > — oo as P^ l LM {x) ~ 1 + (3 c xe~ in full agreement with the analysis of [9]. The 
opposite tail for x — ► oo has a generalized Gumbel-like shape P^ LM {x) oc exp j^p — 2e _ 2~|. 
Note a certain similarity of these two asymptotes to those of the probability density for the 
magnetisation in the low-temperature phase of XY model |21j . 

Moreover, for all temperatures below the transition (3 > (3 C the value of the leading term 
in mi and the shape of the wavefront profile should be frozen to the critical point values, i.e. 
those for (3 = j3 c . Thus, to the leading order mi,(J3 > f3 c ) = c{(3 c )L = m|J§|, whereas the profile 
gpiy) is given by Eq. (|32p for any (3 > (3 C . In the same way as in REM case this fact allows 
one to extract the moments of the partition function everywhere in the low temperature phase 
when G{x) = gp c {x + m* L ). Employing now the critical profile shape Eq. ([32|) and substituting 
x = p Inp we recover the Laplace transform Gp{p) of the probability density of the partition 
function below the transition: 

dZ p Vp >Pc {Zp) e~ pZ " = G p>Pc (p) = Vbp*K x (26p?) , b = e^ m i; j = ^ < 1 . (33) 

This gives us the possibility to calculate negative moments of the partition function as 

_ 2 

( ^ > = r5or* prJG ** W -WM r ( l + ^) r (^)" " >0 (34) 

§ As 0(1) terms in mj, above the transition diverge logarithmically when /3 — > /3 C it is natural to 
expect that at the transition point they should be replaced with constlnL. Actually, the analysis of 
[9] predicts the precise value const = 1/2 at the transition point. Unfortunately, verification of this 
interesting prediction goes beyond the precision of our analysis. 
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where we have used the identity [ 

^^ M< , p = 2'-'< l --'r(l±f^)r(l±^) . (35) 

Substituting here the explicit values of 7 and b, and changing v — ► —v we finally get: 

Jvm* T __}__r 2 ( 1 ^ 



i w^) T v-r)- (36) 

Although we used v < in the course of derivation, a slight modification of the above 
procedure, see [3J, shows that the above expression is valid in a wider region, as long as 
v < c /0 < 1. 

The mean value of the free energy F in the low-temperature phase is found in a similar way 
and the leading order term is simply (F) = —0~ 1 (lnZ) = —m* L . Introducing the probability 
density Vn(f) of / = F + m* L we can now rewrite Eq. (|36p as 

£e*'W>* = FaTI) r2 (! + £). Re s > -/3 (37) 

In particular, similarly to REM case after an analytic continuation s — > is the probability 
density of the free energy for the circular logarithmic model (CLM) can be extracted for any 
> C by inverting the corresponding Fourier transform. The corresponding formula takes a 
form of an infinite series: 

1 f°° 1 ^,9 I \ is 



~df 



^(fi-ZL^f&iff-^P*- (38) 



1 + E- —f ^ /3c/ + --2^(n + l) + ^(l-n§ 



(39) 



where ip{x) = T'(x)/T(x). Exploiting the series expansion for the Macdonald function, see 
e.g. p. 909 of [22J, it is easy to check that in the zero temperature limit — > 00 the free energy 
distribution indeed reduces to the extreme value probability density of the form Eq. ([2]) in full 
agreement with the general relation ([8]). Eq. (|39p shows that the same non-Gumbel behaviour 
holds for the far-left tail / — > —00 of the free energy distribution at any temperature below 
the transition. 



3.3. Conclusion, Discussions and Open Problems 

In the present paper we attempted to investigate some implications of the CLD freezing 
scenario [9] for a particular type of REM-like model with logarithmically correlated random 
potential on a circle. The chosen model seems to be especially attractive due to relatively 
simple expressions for the integer moments of the partition function in the high-temperature 
phase, given by the well-known Dyson Coulomb gas integral. We argue that in such a case 
the Laplace transform of the probability density of the partition function can be efficiently 
recovered. When combined with the freezing scenario this knowledge allows us to continue 
the Laplace transform to the low-temperature phase. We first check that the method indeed 
works for the standard REM example by recovering the well-known, yet nontrivial Gardner- 
Derrida formulae [3] for the moments of the partition function below the freezing point. The 
same method is then applied to the logarithmic model in question. In particular, we are 
able to recover the full distribution of the lowest minimum in the potential, Eq.([2]), and this 
extreme-value statistics is manifestly non-Gumbel. 
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Although we think our results are supported by rather convincing arguments, the 
calculations are very essentially based on a few plausible but not yet fully verified assumptions. 
As such, mathematically our conclusions have the status of well-grounded conjectures. It would 
be certainly very desirable to find alternative ways of investigating the model, as well as to 
perform accurate numerical verification of the precise form of the extreme value statistics. 
Another open problem is the universality of our result, Eq. ([2]), for logarithmically correlated 
random variables, in particular the shape of the right tail (see [17]). We hope our results 
provide enough incentive for further research in this direction. 

Finally, it might be useful to provide an alternative view on our choice of 
the logarithmically correlated potential, Eq. (|15p . By employing the known identity: 
— In (4 sin 2 Xl ^ X2 ) = 2 J2kLi 7 cos K x i ~ x 2) we see that the covariance function (fl~5j) represents, 
in fact, a 27r-periodic real-valued Gaussian random process V(x) = Ya^=i (yi edx + e~ dx ^j 

with a self-similar spectrum (viv m ) = g 2 l~^ 2H+1 ^5i m characterised by the particular choice of 
the Hurst exponent H = 0. Such a process therefore represents a version of the so-called 1/f 
noise. To this end it is worth mentioning that the extreme value statistics of the "roughness" 
associated with 1/f noise was investigated in [23J, and found to be of Gumbel form. 
Acknowledgements. YF acknowledges support of this work by the Leverhulme Research 
Fellowship project "A single particle in random energy landscapes". Discussions with M. 
Feigelman, S. Nechaev, V. Vargas and V. Yudson are acknowledged with thanks. 

Appendix 

The positive integer moments of the distribution V(Z), see Eqs. (f24"|) - (f26|) . are given in the 
high-temperature phase f3g < 1 by the sum of two terms: 

(Z n ) = m { < ] + m { > ] . (40) 

The first contribution corresponding to Eq. (|25p is given by 

m (n) = ^* V< ( Z)Z n dz = Z n j°° r -/?Vn e -r dr ; B=(jZ\ W? . (41) 

_____ 1 m Q^q^\\^\M 

In the limit InM — > oo we have from (|24|) and (|25[) B oc e P 2 9 2 — ► in view of 

(3 2 g 2 < 1. After a simple calculation we find 

m < =\ i z ^ z n ~w^ n> i (42) 

As to the second contribution, a saddle-point analysis justified by InM 3> 1 shows that : 

roo i f(l) (l+2n--_0 lnM , 

m (n) = e n^Z V>{Z) dZ w 2V7r lnM(l-^^) ^ J ' " < (43) 

Jz. [ /(nW)e mM(i+^„») > n> _i_ j 

(n) (n) 

Comparing the two contributions m> and within the high-temperature phase 



(i) m<^ 3> m^ as long as 1 < n < 737-7. Indeed 



(3g < 1 we see that 

(n) „ (n) 

n(l + (3 2 g 2 ) -{l + 2n- J-j) = (1 - p 2 g 2 ) (-jL - n) > 

which implies 



» 
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(ii) If n > -g^p- we have m< <C m> , as in this case 

4"» ~ zf^Z^ ~ e ( 1+2n w) lnA/ « ro W ~ e (i+W)inM i 
which follows from 

( lt «» ! )-(l + 2»-^) = (j 9 »-i) ! >». 

Accounting for Eq. (|23|) and the definition of Z e in Eq. (j25j) we indeed see that the moments 
{Z n ) coincide for InM > 1 with the expressions for the partition function moments (|2ip . 
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